Many-body effects in TiSe2: can GW describe an excitonic insulator? 
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We present both theoretical ab initio GW and experimental angle-resolved photoemission 
(ARPES) and scanning tunneling (STS) spectroscopy results on TiSe2. With respect to the density- 
functional Kohn-Sham metallic picture, the many-body GW self-energy leads to a « 0.2 eV band 
gap insulator consistent with our STS spectra at 5 K. The band shape is strongly renormalized, with 
the top-of-valence moved towards a circle of points away from P, arising in a mexican hat feature 
typical of an excitonic insulator. Our calculations are in good agreement with experiment. 

PACS numbers: 71.15.Qc, 71.20.-b, 71.30.+h, 71.35.-y 
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Introduction Crystalline solids are classified as either 
insulators or metals, depending on the existence or ab- 
sence of a finite energy gap in the electronic excita- 
tion spectrum. Within the independent particle picture, 
metallic or insulating behavior is explained in terms of 
partially or completely filled valence bands, but elec- 
tronic correlations may overcome this ordinary paradigm 
and lead to new mechanisms for transforming a metal 
into an insulator or viceversa. From the early works of 
Mott, Kohn, and Hubbard [iHl], several systems show- 
ing such exotic behavior have been supposed to exist in 
nature. One of them is the so-called excitonic insula- 
tor [H, 0-0] , first predicted and studied in the works of 
Mott. Kopaev, Keldysh and Kohn and not to be con- 
fused with the more celebrated Mott-Hubbard insulator. 
Here we study the possibility of an excitonic insulator in 
titanium diselenide (TiSe2) 0], focusing on the question 
whether ab initio many-body GW theory [^, can de- 
scribe photoemission spectra of such an excitonic phase. 

At a non-interacting one-electron level, an excitonic in- 
sulator would present a semimetallic or semiconducting 
band structure with sufficiently small (typically few tens 
of meV) overlap or gap, and reduced number of free car- 
riers. Consequently, the Coulomb interaction between 
particles is only weakly screened. Electrons and holes 
can then spontaneously bind into non-conducting exci- 
tons and form a new ground state of lower energy than 
the normal phase. Exciton condensation may also lead to 
the formation of charge-density waves (CDW) of purely 
electronic origin. On the semimetal side, model descrip- 
tions of the system are based on a BCS-like approach 
where the role of electron-electron Cooper pairs is taken 
by electron-hole excitons. Similar to BCS superconduc- 
tors, an energy gap opens of the same order of magnitude 
as the binding energy of the pair [llj (depending on the 
underlying band structure, e.g. effective mass anisotropy 
and mult i valley effects). 



So far there is no clear-cut experimental proof of the 
excitonic insulator existence in nature. Nevertheless, sev- 
eral experimental obervations poin t to TiSe2 as one of 
the most probable candidates 12l-ll7l|. Indeed a CDW 
has been observed in the ground state of TiSe2, like in 
other transition-metal dichalcogenides, but, so far, there 
is no consensus on the interpretation of the ground state 
and its CDW. In particular, the exact role played by ex- 
citons is still controversial. Either a Jahn- Teller effect 



19| or a correlation mechanism leading to the exci- 
tonic insulator 13|, [ij] have been invoked. Recently, it 
was suggested that the ground state and the resulting 
periodic lattice distortion (PLD) that occurs below 202 
K are primary due to electron-phonon coupling, which 
is only enhanced by the presence of incoherent excitons 
[20t . On the other hand, photoemission spectra are con- 
sistent with the excitonic insulator scenario |15| , and the 
interaction between an exciton condensate and phonons 
can reproduce the observed atomic displacements in the 
PLD with good agreement [l7| . 

In this work we provide theoretical evidence from ab 
initio many-body GW calculations supporting the exci- 
tonic insulator scenario in TiSe2 and verify our calcula- 
tions against angle-resolved photoemission spectroscopy 
(ARPES) measurements and scanning tunneling spectra 
(STS). STS at 5 K clearly show a bandgap of ~ 0.15 
eV and point to TiSe2 as a semiconductor. With re- 
spect to the one-electron mean-field density-functional 
theory (DFT) band structure, where TiSe2 is a metal, the 
GW many-body self-energy opens a band gap of ~ 0.2 
eV, thus leading to an insulator. Furthermore, the band 
shape at high symmetry points is strongly renormalized. 
In particular, the top-of-valence is moved from P toward 
a circle of points away from P, building up a mexican hat 
feature characteristic of an excitonic insulator, as pre- 
dicted by Kohn [1] . ARPES spectra indicate a flattening 
of the topmost valence bands before P, compatible with 
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FIG. 1: Left panel: STS dl/dV experimental spectra (green 
circles) vs GW (shifted to account for finite impurity dop- 
ing of the sample, black curve) and DFT (zeroed at Ef, 
red curve) DOS for TiSe2. Right panel: Energy- loss func- 
tion. Green dots: experiment [23| (derived from ellipsome- 
try); black curve: present theoretical RPA result. 



the mexican hat picture. The GW bands we obtain are 
in good agreement with ARPES measurements, and also 
our STS spectra compare more favorably with the GW 
density-of-states (DOS) than with DFT. The calculated 
dielectric function describing the screening properties of 
the system, is in good agreement with previously mea- 
sured energy- loss spectra (EELS) The overall pic- 
ture arising from our results indicates a possible ground 
state instability in TiSe2 of electronic origin. Further, 
our calculations show that bare GW is able to repro- 
duce the photoemission band structure of systems with 
strong excitonic renormalizations, although, in general, 
vertex corrections beyond GW are necessary to properly 
describe excitons in optical spectra [131 ■ 

Theory Our first step is a standard, numerically well 
converged [IB] ground-state DFT-LDA calculation of 
the total energy, the electronic density, and the lat- 
tice parameters within the non-reconstructed primitive- 
hexagonal crystal structure. The DFT calculation pro- 
vides the Kohn-Sham (KS) electronic structure that we 
use for calculating the Green's function G and the dy- 
namically screened interaction W{uj) = e~^{uj)v — de- 
fined as the bare Coulomb interaction v screened by the 
dynamical dielectric function s~^{uj) — , both entering 
into the self-energy in the G W approximation [13] , 
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Within a perturbative approach to first-order in S — 
w^J^*, we then calculate the GW quasiparticle energies, 
expanding S around the KS energies e^^ , taken as a zero- 
th order approximation for the true quasiparticle ener- 
gies. 
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where Z = (1 - dY.^"^ /duj\ _ ks) 



In the standard GW approach, the dielectric function 
e~^{u!) is calculated in the RPA approximation. Since 
screening is of crucial importance for exciton formation, 
we have checked the RPA energy-loss function — 3e~^(w) 
against the results of ellipsometry experiments (see 
Fig. [1]). From the good agreement with experiment, we 
are confident that the used RPA screening accurately de- 
scribes the physical situation. Notice that the effective 
screening of our ab initio RPA calculation is much weaker 
compared to the Lindhard function, used in the model 
calculation of Ref. [l^, so that we expect that exciton 
formation is even more favored. 

Experiment Photon-energy dependent ARPES mea- 
surements were carried out at the SIS beamline of the 
Swiss Light Source synchrotron, using a Scienta SES- 
2002 spectrometer with an overall energy resolution bet- 
ter than 10 meV, and an angular resolution better than 
0.5 deg. The data shown here have been collected at 65 
K on two samples from the same batch with a slight Ti 
overdoping, resulting in a small population of the elec- 
tron pocket at the L point and a consequent shift of the 
valence band. Scanning tunneling spectroscopy (STS) 
spectra were measured at low temperature by use of a 
OMICRON LT-STM with a lock-in technique. The esti- 
mated energy resolution is about 10 meV. 

Results: STS and GW band gap The results of our 
STS measurements at 5 K are given in Fig. [TJ We plot 
the differential conductance (dl/dV) as a function of the 
bias V. For small enough biases and slowly varying tun- 
neling matrix elements, the differential conductance gives 
an image of the local density-of-states (DOS) that can 
be compared with calculations. At small bias the ex- 
perimental dl/dV vanishes, indicating a gap of ~ 0.15 
eV. With respect to the long debated TiSe2 semimetal- 
semiconductor question 2M j our STS spectra, the first to 
date to our knowledge, provide a clear indication of semi- 
conducting character, at least at low temperature. The 
STS gap is in quantitative agreement with recent ARPES 
experiments [l3, 14, l^, 21, 22|. In Fig. [Done can see that 
the calculated DFT DOS never vanishes, thus indicat- 
ing a metal. On the other hand, the GW DOS presents 
a well defined gap, like in STS, with just only a slight 
overestimation (~ 0.2 eV). Contrary to DFT, GW found 
TiSe2 a very small gap insulator, in agreement with the 
experiment. The left-valence and right-conduction DOS 
profiles are also well reproduced by GW. 

Results: the band plot In Fig. [2]a we report the band 
plot of TiSe2 based on DFT-LDA and GW calculations, 
together with experimental ARPES data. According to 
the DFT KS electronic structure, TiSe2 is metallic with 
a band overlap of 0.4 eV between the top of the pocket 
of holes at F and the bottom of the pocket of electrons 
at L. The GW electronic structure then leads to a gen- 
eral shift of the bands, negative for the lowest and pos- 
itive for the highest, and opens a band gap of around 
0.2 eV. Therefore, many-body effects, as accounted by 
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FIG. 2: TiSe2 band plot (a panel) and its zoom (d) at F (M-F-A directions). Red lines: DFT-LDA KS electronic structure; 
yellow dashed lines; DFT plus spin-orbit correction (SOC); black dots and lines: GW band plot; violet dots: GW results 
including Ti 3s3p electrons in valence; blue diamonds and green squares: sample 1 and 2 ARPES experimental data at T = 65 
K. Panel b: ARPES signal I{k,E); and c: its derivative (fl/dE^ taken at T = 65 K around F, both in the F-M direction. 



the GW approximation, turn the KS metal into a small 
gap insulator. As this happens also in other systems, 
e.g. germanium, it is not a novelty in itself. However, 
here the important modification of the band shape shown 
in Fig. is quite unusual: The top-of-valence is moved 
away from F in the plane perpendicular to towards a 
circle of points around F. Similar reshaping also appears 
for the lowest two empty bands at F and at A, and, less 
pronounced, for two lower bands at A and M. Such mexi- 
can hat features have been predicted by Kohn Q to occur 
in excitonic insulators. Up to our knowledge, our results 
are the first observation of this characteristic shape from 
ab initio G W calculations, pointing to an excitonic origin 
of the insulating nature of TiSe2. 

Our (Fig. (2?*) and previous ARPES experiments [13, 
21[ show a weakening of the signal in the neighbourhood 
of F, so that the band structure cannot be clearly re- 
solved there. However Kidd et al. first remarked a loss 
of parabolic shape of the Se 4p band, together with a flat- 
tening (see Fig. lb and text in Ref. [14|) occurring at fc ~ 
0.05/FK which is confirmed also by our ARPES spectra 
(Fig. Our GW bands (Fig. ^) present a loss of 

parabolic shape and a maximum, i.e. a zero of the first 
derivative de/dk, and so a flattening at fc ~ 0.05/FM, 
perfectly compatible with Kidd et al. and our ARPES 
findings. A plot of the second derivative d'^I/dE^ of the 
ARPES signal (Fig. ^) presents further anomalies: at F 
we can observe a peak at ~ —0.5 eV, to be interpreted 
as the light electron F^ band; and then two peaks at 
~ —0.3 and ~ —0.2 eV (red squares in Fig. [21s), proba- 
bly corresponding to a split of the Fg" Sc 4p degenerate 



bands at F. Since those bands top at i? ~ —0.2 eV and 
at £^ ~ —0.15 eV (red circles in Fig. [SJc), we infer that a 
band bending should occur in order to achieve the low- 
est energies at F. This points to a possible experimental 
confirmation of the GW mexican hat scenario. 

For the rest of the band plot (Fig. ^jp, and d), exper- 
imental ARPES points agree much better with the GW 
bandplot than with DFT KS, both in the FM and in FA 
directions. DFT presents empty Se 4p bands around F 
and along FA, giving rise to a pocket of holes at F. This 
is in striking contrast with ours (Fig.[2j)) and all previous 
ARPES results [US El, [13 which found occupied Se 
4p bands and no pocket of holes, exactly like in our GW 
calculation. Although the question was long debated, the 
most recent ARPES experiments [H, [13, [11, [HI, [llj seem 
to indicate an indirect band gap of ~ 150 meV between 
the Se 4p bands at F and the "emerging" (above Ep and 
just only termally occupied) Ti 3d-derived band at L. 
This scenario is in good agreement with our GW band 
structure and in contrast with DFT and oldest ARPES 
experiments which rather found a band overlap. Resid- 
ual mismatches of our GW and ARPES bands, like the 
split of the two Se Ap bands along FA, can be solved 
by including the spin-orbit coupling (SOC), as one can 
see by comparing DFT and DFT SOC bands in Fig. 
SOC and relativistic corrections have also important ef- 
fects in band-crossing avoidings. Finally, the mismatch 
in the position of the flat Ti 3d band at -2.5 eV is due 
to the lack of scmicore electrons. Inclusion of semicore 
Ti 3.s3p electrons into the pseudopotential valence shell 
(violet dots in Fig. ^jp,) improves the position of the 3c? 
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FIG. 3: Exchange {nk\'Ex\nk) (black circles) and correlation 
{nk\T,^'^ \nk) (red diamonds) real part contributions to the 
self-energy matrix elements, for the 3 topmost occupied (filled 
symbols) and 3 lowest empty (empty symbols) bands along 

M-r. 



bands with respect to the experiment. 

Microscopic analysis of the Mexican hat scenario In 
Fig. 13] where we report the separate contributions of the 
diagonal matrix elements to the correlation {nk\Y]2^ \nk) 
and exchange {nk\Yix\nk) self-energy, entering into the 
GW correction to the KS energies, for the 3 topmost 
occupied and the 3 lowmost empty states at F. It is evi- 
dent that the mexican hat renormalization is dominated 
by exchange effects. Bare exchange introduce changes of 
2 eV between F and the new top-of-valence (TOV) circle 
of maxima, while screening effects included in the corre- 
lation part reduce this effect to a total value of around 
0.3 eV. Similar conclusions hold also for the two low- 
most empty states at F, although with a softer renormal- 
ization. If the mexican hat many-body renormalization 
that we see in TiSe2 is really associated to its excitonic 
insulator nature, then it is already captured at a low 
level of approximation on the self-energy and there is no 
need to improve on correlation. In general, one expects 
that vertex corrections beyond the GW approximation 
are required to properly describe excitons, e.g. within 
the Bethe-Salpeter equation approach [2J] to describe the 
spectrum of neutral excitations as in optical absorption. 
However, in contrast to optical absorption, ARPES mea- 
sures charged excitations, and it seems that a lower level 
approximation, as the GW used here, can describe ex- 
citonic effects in ARPES spectra. A Feynman diagram 
expansion of the GW Green's function would show that 
hole (electron) channels are taken into account by the 
GW renormalization of the electron (hole) propagator. 

In order to understand the mechanism and to identify 
the states responsible of the mexican hat effect, we stud- 
ied the matrix elements of the electric dipole operator 
D = e"*''", directly entering into the exchange operator: 



4^ ^ „ KV'n'fe'|e-''('-'')lV^nfc)|' 



(1) 



with f„k Fcrmi-Dirac occupation numbers. We studied 
in particular the difference of matrix elements calculated 



at A: = F and those at k = TOV (the top of the mex- 
ican hat), as a function of n' and k' . We found that 
the mexican hat feature for the topmost n = Se 4p3/2 
band is mostly conjured by exchange with states on bands 
n' = Se 4p3/2 and 4pi/2j and at k' close to the Brillouin 
zone center and rather in the k^ky plane (we checked 
the FM direction). On the other hand we have found 
negligible the contribution coming from exchange with 
states at k' along the kz direction (FA direction), and 
also from states faraway from A: ~ F where the mexican 
hat feature occurs — the Brillouin zone boundary for 
instance — . The latter are suppressed by the \k' — k\^ 
factor at the denominator in Eq. ([T]). Our analysis seems 
also to exclude that the mexican hat is an artifact of 
our GqWq non self-consistent calculation. Indeed the Ti 
3d-derived band, which is partially filled in DFT and be- 
comes empty in GW, plays no role in the renormalization. 
In the next iterations the increased population of the Se 
Ap bands, mostly responsible of the renormalization, and 
the opening of a bandgap, resulting in a reduction of the 
electron-hole screening, should lightly further strengthen 
the renormalization, or at most not change the picture. 
Overall, our analysis qualitatively confirms the physics 
emerging from the model calculation of Ref. up to the 
possibility of a CDW. 

Finally, our calculation shows that the mexican hat 
effect is already present for the higher symmetry non- 
distorted atomic structure and/or without invoking the 
role of ionic degrees of freedom such as those involved 
in the electron-phonon interaction. If the mexican hat 
is directly associated to the excitonic insulator nature of 
the system, it will indicate an instability of the ground 
state of purely electronic origin, without the necessity of a 
phonon-driven ionic mechanism. Hence, the experimen- 
tally observed 2x2x2 distorted reconstruction [I3 should 
rather be regarded as a consequence, following the exci- 
tons condensation and the formation of a CDW, rather 
than a cause. 

Conclusions We presented both theoretical GW and 
new experimental ARPES and STS results on TiSe2 . The 
many-body GW self-energy opens a band gap of w 0.2 
eV, and strongly renormalize the band shape with the 
top-of-valence moved toward a circle of points away from 
F, arising in a mexican hat feature typical of an excitonic 
insulator. The calculations are in good agreement with 
our ARPES and STS data, which are compatible with 
the mexican hat scenario, as well as EELS experiments. 
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